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ABSTRACT 

This paper deals with the comparison of planar parallel ma- 
nipulator architectures based on a multi-objective design opti- 
mization approach. The manipulator architectures are compared 
with regard to their mass in motion and their regular workspace 
size, i.e., the objective functions. The optimization problem is 
subject to constraints on the manipulator dexterity and stiffness. 
For a given external wrench, the displacements of the moving 
platform have to be smaller than given values throughout the 
obtained maximum regular dexterous workspace. The contri- 
butions of the paper are highlighted with the study of 3-PRR, 
3- RPR and 3-RRR planar parallel manipulator architectures, 
which are compared by means of their Pareto frontiers obtained 
with a genetic algorithm. 



INTRODUCTION 

The design of parallel kinematics machines is a complex 
subject. The fundamental problem is that their performance 
heavily depends on their geometry UJ and the mutual depen- 
dency of almost all the performance measures. This makes the 
problem computationally complex and yields the traditional so- 
lution approaches inefficient. As reported in 1|2J, since the perfor- 
mance of a parallel manipulator depends on its dimensions, the 
latter depend on the manipulator application(s). Furthermore, 
numerous design aspects contribute to the Parallel Kinematics 
Machine (PKM) performance and an efficient design will be one 
that takes into account all or most of these design aspects. This is 



an iterative process and an efficient design requires a lot of com- 
putational efforts and capabilities for mapping design parame- 
ters into design criteria, and hence turning out with a multiobjec- 
tive design optimization problem. Indeed, the optimal geomet- 
ric parameters of a PKM can be determined by means of a the 
resolution of a multiobjective optimization problem. The solu- 
tions of such a problem are non-dominated solutions, also called 
Pareto-optimal solutions. Therefore, design optimization of par- 
allel mechanisms is a key issue for their development. 

Several researchers have focused on the optimization prob- 
lem of parallel mechanisms the last few years. They have come 
up either with mono- or multi-objective design optimization 
problems. For instance, Lou et al. Hill) presented a general ap- 
proach for the optimal design of parallel manipulators to maxi- 
mize the volume of an effective regular-shaped workspace while 
subject to constraints on their dexterity. Hay and Snyman IT] 
considered the optimal design of parallel manipulators to obtain 
a prescribed workspace, whereas Ottaviano and Ceccarelli |l5]|6l 
proposed a formulation for the optimum design of 3-Degree- 
Of-Freedom (DOF) spatial parallel manipulators for given po- 
sition and orientation workspaces. They based their study on the 
static analysis and the singularity loci of a manipulator in order 
to optimize the geometric design of the Tsai manipulator for a 
given free-singularity workspace. Hao and Merlet Q discussed 
a multi-criterion optimal design methodology based on interval 
analysis to determine the possible geometric parameters satisfy- 
ing two compulsory requirements of the workspace and accuracy. 
Similarly, Ceccarelli et al. jH) dealt with the multi-criterion op- 
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timum design of both parallel and serial manipulators with the 
focus on the workspace aspects, singularity and stiffness prop- 
erties. Gosselin and Angeles ||9][T0l analyzed the design of a 
3-DOF planar and a 3-DOF spherical parallel manipulators by 
maximizing their workspace volume while paying attention to 
their dexterity. Pham and Chen ifTTl suggested maximizing the 
workspace of a parallel flexible mechanism with the constraints 
on a global and uniformity measure of manipulability. Stamper et 
al. lfT2l used the global conditioning index based on the integral 
of the inverse condition number of the kinematic Jacobian matrix 
over the workspace in order to optimize a spatial 3-DOF trans- 
lational parallel manipulator. Stock and Miller [13| formulated 
a weighted sum multi-criterion optimization problem with ma- 
nipulability and workspace as two objective functions. Menon 
et al. IIT4II used the maximization of the first natural frequency 
as an objective function for the geometrical optimization of the 
parallel mechanisms. Similarly, Li et al. flSl proposed dynam- 
ics and elastodynamics optimization of a 2-DOF planar parallel 
robot to improve the dynamic accuracy of the mechanism. They 
proposed a dynamic index to identify the range of natural fre- 
quency with different configurations. Krefft |fT6l also formulated 
a multi-criterion elastodynamic optimization problem for paral- 
lel mechanisms while considering workspace, velocity transmis- 
sion, inertia, stiffness and the first natural frequency as optimiza- 
tion objectives. Chablat and Wenger [ 17] proposed an analytical 
approach for the architectural optimization of a 3-DOF transla- 
tional parallel mechanism, named Orthoglide 3-axis, based on 
prescribed kinetostatic performance to be satisfied in a given 
Cartesian workspace. 

Most of the foregoing research works aimed to improve the 
performance of a given manipulator and the comparison of vari- 
ous architectures for a given application or performance has not 
been considered. In this paper, the mechanisms performance are 
improved over a regular shaped workspace for given specifica- 
tions. As a result, we propose a methodology to deal with the 
multiobjective design optimization of PKMs. The size of the 
regular shaped workspace and the mass in motion of the mech- 
anism are the objective functions of the optimization problem. 
Its constraints are determined based on the mechanism accuracy, 
assembly and the conditioning number of its kinematic Jacobian 
matrix. The proposed approach is applied to the optimal design 
of Planar Parallel Manipulators (PPMs) with the same mobility 
and set of design parameters. The non-dominated solutions, also 
called Pareto-optimal solutions, are obtained by means of a ge- 
netic algorithm for the three architectures and finally a compari- 
son is made between them. 



MANIPULATORS UNDER STUDY 

Figure |l(a)| -(c) illustrate the architectures of the planar par- 
allel manipulators (PPMs) under study, which are named 3-PRR, 
3-RPR and 3-RRR PPMs, respectively. Other families of PPMs 



are described in IJl- Here and throughout this paper, R, P, R 
and P denote revolute, prismatic, actuated revolute and actuated 
prismatic joints, respectively. The manipulators under study are 
composed of a base and a moving platform (MP) connected by 
means of three legs. Points Ai,A2 and A3, (Ci, C2 and C3, respec- 
tively) lie at the corners of a triangle, of which point O (point P, 
resp.) is the circumcenter. Each leg of the 3-PRR PPM is com- 
posed of a P, a R and a R joint in sequence. Each leg of the 
3-RPR PPM is composed of a R, a P and a R joint in sequence. 
Likewise, each leg of the 3-RRR PPM is composed of three R 
joints in sequence. The three P joints of the 3-PRR and the 3- 
RPR PPMs are actuated while the first R joint of each leg of the 
3-RRR PPM is actuated. 

^6 and -J^p are the base and the moving platform frames 
of the manipulator. In the scope of this paper, and are 
supposed to be orthogonal. is defined with the orthogonal 
dihedron {Ox,Oy), point O being its center and Ox parallel to 
segment A 1A2. Likewise, J^p is defined with the orthogonal di- 
hedron {PX,PY), point C being its center and PX parallel to seg- 
ment C1C2. The manipulator MP pose, i.e., its position and its 
orientation, is determined by means of the Cartesian coordinates 
vector p = [px,Py]^ of operation point P expressed in frame ^i, 
and angle 0, namely, the angle between frames and J^p. 

The geometric parameters of the manipulators are defined 
as follows: (i) R is the circumradius of triangle A1A2A3 of cir- 
cumcenter O, i.e., R ~ OAj; (ii) r is the circumradius of triangle 
C1C2C3 of circumcenter P, i.e., r ~ PCt, i = 1, ... ,3; (iii) Lf, is 
the length of the intermediate links, i.e., Li, ~ BiCt for the 3- 
PRR PPM. Lj, is also the maximum displacement of the pris- 
matic joints of the 3-RPR PPM. Similarly, Ly is the length of the 
two intermediate links of the 3-RRR PPM, i.e.. Lb =AiBi = Bid; 
(iv) rj is the cross-section radius of the intermediate links; (v) rp-. 
the cross-section radius of links of the moving platform, the latter 
being composed of three Unks. 

Stiffness Modeling 

The stiffness models of the three manipulators under study 
are obtained by means of the refined lumped mass modeling de- 
scribed in |fT8|. Figures ID to |4] illustrate the flexible models of 
the legs of the 3-PRR, 3-RPR and 3-RRR PPMs, respectively. 
The actuator control loop compliance is described with a 1-dof 
virtual spring and the mechanical compliance of each link with 
a 6-dof virtual spring in each flexible model denoted 0,. Be- 
sides, the moving platform of the manipulators is supposed to be 
composed of three links of length r connected to its geometric 
center P. 

From Fig. [2] the flexible model of the legs of the 3- 
PRR PPM contains sequentially: (i) a rigid link between the 
manipulator base and the actuated joint (part of the base plat- 
form) described by the constant homogeneous transformation 
matrix T^^^,^; (ii) a 1-dof actuated joint, defined by the homo- 



2 



Copyright (c) 2010 by ASME 




geneous matrix function Va(q'Q) where q'^ is the actuated coor- 
dinate; (iii) a 1-dof virtual spring describing the actuator me- 
chanical stiffness, which is defined by the homogeneous ma- 
trix function Ysi {%) where 6q is the virtual spring coordinate 



corresponding to the translational spring; (iv) a 1-dof passive R- 
joint at the beginning of the leg allowing one rotation angle q'2, 
which is described by the homogeneous matrix function V^i {q'j)', 
(v) a rigid leg of length L Unking the foot and the movable plat- 
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FIGURE 4. FLEXIBLE MODEL OF THE 3-RRR PPM'S KINE- 
MATIC CHAINS 

form, which is described by the constant homogeneous transfor- 
mation matrix T^; (vi) a 6-dof virtual spring describing the leg 
stiffness, which is defined by the homogeneous matrix function 
V,2 {9{ ■ ■ ■ 0^), with e{, 9i,, 9l and 0^, 0^ being the virtual 
spring coordinates corresponding to the spring translational and 
rotational deflections; (vii) a 1-dof passive /? -joint between the 
leg and the platform, allowing one rotation angle ^3, which is 
described by the homogeneous matrix function Vri (^3 ) ; (viii) a 
rigid link of length r from the manipulator leg to the geometric 
center of the mobile platform, which is described by the con- 
stant homogeneous transformation matrix Tj.; (ix) a 6-dof vir- 
tual spring describing the stiffness of the moving platform, which 
is defined by the homogeneous matrix function V53 (0^ • • • 6[2), 
0j, 0g, 0g and 0Jq, 0Jj, 0J2 being the virtual spring coordinates 
corresponding to translational and rotational deflections of link 
C,P; (x) a homogeneous transformation matrix T^^^ that char- 
acterizes the rotation from the 6-dof spring associated with link 
CiP and the manipulator base frame. 

As a result, the mathematical expression defining the end- 
effector location subject to variations in all above defined coor- 
dinates of a single kinematic chain ; of the 3-PRR PPM takes the 
form; 



TWVi (qh) V,i (0^) V,i {q\) T[V.2 (0; • • • 0^) 

V.2(^yTtV.3(0f ••0l2)n„, 



(1) 



Similarly, the mathematical expressions associated with the 
kinematic chains of the 3-RPR and 3-RRR PPMs are obtained. 

From HH, the kinetostatic model of the if/j leg of the X- 
PPMs can be reduced to a system of two matrix equations, 
namely. 



(2) 
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the virtual springs and is the one related to the passive joints. 
' describes the compliance of the virtual springs. 
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where K'^^., is the 1x1 stiffness matrix of the itii actuator, KJ^-^^^, is 
the 6x6 stiffness matrix of the intermediate link for the 3-PRR 
and 3-RPR PPMs while Kj.,,^.^ and K;,.„^^ are the 6 x 6 stiffness 
matrices of the first and second intermediate links of the ith leg 
of 3-RRR PPM. K^^^ is the 6 x 6 stiffness matrix of the ith Unk 
of the moving platform. The compliance matrix of each link is 
expressed by means of the stiffness model of a cantilever beam, 
namely. 
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(4) 



L being the length of the corresponding link, A is its the cross- 
sectional area, i.e., A = nrj for the links of the manipulators legs 
and A = Ttr^ for the links of the moving platform. ly and 4 are the 
polar moments of inertia abouty and z axes, resp. ly = Iz = ?rrj/4 
for the links of the manipulators legs and ly — Iz = for the 

links of the moving platform. I^=I-+Iy is the polar moment of 
inertia about the longitudinal axis of the link. E and G are the 
Young and shear moduli of the material. 

Accordingly, the Cartesian stiffness matrix K, of the ith leg 
defining the motion-to-force mapping is obtained from Eq. (|2]i. 



where X stands for 3-PRR, 3-RPR or 3-RRR. The sub-matrix 

— 1 ■ T 

^'e\x ^ Je|x^e|x J'eix describes the spring compliance rela- 
tive to the geometric center of the moving platform, and the sub- 
matrix takes into account the passive joint influence on the 
moving platform motions. Jg is the Jacobian matrix related to 



f,=K,5t; 



(5) 



with f; being the wrench exerted on the ith leg of the manipulator 
and at the geometric center of the moving platform while 5t; is 
the small-displacement screw of the moving-platform. 
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Finally, the Cartesian stiffness matrix K of the manipulator 
is found with a simple addition of the three K, matrices, namely, 

3 

K = £K,- (6) 

!=1 

MULTIOBJECTIVE OPTIMIZATION PROBLEM 

A multiobjective optimization problem (MOOP) is formu- 
lated in this section in order to compare 3-PRR, 3-RPR and 
3-RRR PPMs. In scope of this study, the manipulators are 
compared with regard to their mass in motion and their regular 
workspace size, i.e., the two objective functions of the MOOP, 
defined below. Moreover, the MOOP is subject to constraints on 
the manipulator dexterity and stiffness. It means that for a given 
external wrench, the displacements of the moving platform have 
to be smaller than given values throughout the obtained maxi- 
mum regular dexterous workspace. 

Objective Functions 

Mass in Motion of the Manipulators The compo- 
nents in motion of the manipulators are mainly their moving plat- 
form and the links of their legs. As a consequence, the mass in 
motion for the three PPMs under study is expressed as follows: 

'WPRR = 'imiink + mpf (7a) 
ffiRpR = 3 mii„i, + nipf (7b) 
WRRR = 6 miink + mpf (7c) 

mii„ii is the mass of links of the legs and are supposed to be the 
same while nipj is the mass of the moving platform. The mass of 
the prismatic or revolute actuators does not appear in Eqs. jTaj i- 
(c) as it is supposed to be fixed for the 3-PRR PPM and close to 
the base for the 3-RPR PPM. 

nipf = nPprv (8a) 
miink^nr]Lv (8b) 

where v is the material density. 

Finally, the first objective function of the MOOP is ex- 
pressed as: 

/i (x) = mx min (9) 

X being the vector of design variables, i.e., the geometric parame- 
ters of the manipulator at hand, and X stands for 3-PRR, 3-RPR 
or 3-RRR. 



Regular workspace size The quahty of the manipula- 
tor workspace is of prime importance for the design of Parallel 
Kinematics Machines (PKMs). It is partly characterized by its 
size and shape. Moreover, the lower the amount of singularities 
throughout the workspace, the better the workspace for continu- 
ous trajectory planning. The workspace optimization of parallel 
manipulators can usually be solved by means of two different 
formulations. The first formulation aims to design a manipulator 
whose workspace contains a prescribed workspace and the sec- 
ond one aims to design a manipulator, of which the workspace 
is as large as possible. However, maximizing the manipulator 
workspace may result in a poor design with regard to the manip- 
ulator dexterity and manipulability lfT2l[T9l . This problem can 
be solved by properly defining the constraints of the optimiza- 
tion problem. Here, the multiobjective optimization problem of 
PPMs is based on the formulation of workspace maximization, 
i.e, the determination of the optimum geometric parameters in 
order to maximize a regular-shaped workspace. 

In the scope of the paper, the regular-shaped workspace is 
supposed to be a cylinder of radius for which at each point 
a rotation range = 20 ° of the moving-platform about the Z- 
axis has to be reached. Figure|5]illustrates such a regular-shaped 
workspace, whose Xc, jc and 0^ are its center coordinates and the 
rotation angle of the moving-platform of the manipulator in the 
home posture. 




FIGURE 5. A REGULAR-SHAPED WORKSPACE 



Consequently, in order to maximize the manipulator 
workspace, the second objective of the optimization problem can 
be written as: 

/2(x)=/?,,^max (10) 

Constraints of the Optimization Problem 

The constraints of the optimization problem deals with the 
geometric parameters, the dexterity and the accuracy of the ma- 
nipulators. Moreover, the constraints have to be defined in order 
to obtain a singularity-free regular-shaped workspace. 
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Constraints on the Geometric Parameters For the 

three PPMs under study, the kinematic constraints are handled 
with their inverse kinematics. It means that the inverse kine- 
matics is solved in order for the postures of the PPM to belong 
to the same working mode throughout the manipulator regular- 
shaped workspace. Besides, for the 3-PRR PPM, the lower and 
upper bounds of the prismatic lengths p, are defined such as 
< p, < \/3R in order to avoid collisions. To obtain feasible dis- 
placements of the prismatic joints, the range of the 3-RPR PPM 
is defined such that L/2 < p, < L. 

Constraint on thie IVIanipulator Dexterity The ma- 
nipulator dexterity is defined by the condition number of its kine- 
matic Jacobian matrix. The condition number K'f(M) of a m x n 
matrix M, with ni < «, based on the Frobenius norm is defined 
as follows 

Kf{M) = -A/tr(M^M)tr[(M^M)-i] (11) 

Here, the condition number is computed based on the Frobenius 
norm as the latter produces a condition number that is analytic 
in terms of the posture parameters whereas the 2-norm does not. 
Besides, it is much costlier to compute singular values than to 
compute matrix inverses. 

The terms of the direct Jacobian matrix of the three PPMs 
under study are not homogeneous as they do not have same units. 
Accordingly, its condition number is meaningless. Indeed, its 
singular values cannot be arranged in order as they are of dif- 
ferent nature. However, from ||20l and | |2T) . the Jacobian can be 
normalized by means of a normalizing length. Later on, the con- 
cept of characteristic length was introduced in ll22l in order to 
avoid the random choice of the normalizing length. For instance, 
the previous concept was used in 1231 to analyze the kinetostatic 
performance of manipulators with multiple inverse kinematic so- 
lutions, and therefore to select their best working mode. 

Accordingly, for the design optimization of the three PPMs, 
the minimum of the inverse condition number fc^* (J) of the 
kinematic Jacobian matrix J is supposed to be higher than a pre- 
scribed value, say 0.1, throughout the regular-shaped workspace, 
for any rotation of its moving-platform, i.e., 

mm (k--' (J)) > 0.1 (12) 

Constraints on thie moving-platform pose errors 

The position and orientation errors on the moving-platform are 
evaluated by means of the stiffness models of the manipulators. 
Let {dx, 8y, dz) and (5^.,, 8(j)y, 8(j)z) be the position and orien- 
tation errors of the moving-platform subject to external forces 
{Fx, Fy, F^) and torques (t^, Ty, T,). The constraints on the pose 



errors on the moving-platform are defined as follows: 

dx < Sx'""" 5y < 8y'"'" 5z < 

(13) 

S(px < ^C'"" 5f . < 50™'-'^ 5(1)- < 

(5x"™, Sy"'"", 5z""") being the maximum allowable position 
errors and (50™", 50""", 50""") the maximum allowable ori- 
entation errors of the moving-platform. These accuracy con- 
straints can be expressed in terms of the components of the mech- 
anism stiffness matrix and the wrench applied to the moving- 
platform. Let us assume that the accuracy requirements are: 



y^Sx^- + 5y2 < 0.0001m (14a) 
5z< 0.001m (I4b) 
50, <ldeg (14c) 

If the moving -platform is subject to a wrench whose components 
are \\F-f.y\ \ =F;= lOON and T- = 100 Nm, then the accuracy con- 
straints can be expressed as: 

k^" > \\Fx.y\\ /VSx^TSf = 10^ N.m-' (15a) 
kf" > F:J8z = 10^ N.m"^ (15b) 

k'r > T,/50, = N.m.rad-1 (15c) 

;r/l50 

Design Variables of the Optimization Problem 

Along with the above mentioned geometric parameters (R, r, Lf,) 
of the PPMs, the radius rj of the circular-cross-section of the 
intermediate bars defined and the radius rp of the circular-cross- 
section of the platform bars are considered as design variables, 
also called decision variables. As a remainder, the moving- 
platform is supposed to composed of three circular bars of length 
r. 

As there are three PPMs under study, the PPM type is an- 
other design variable that has to be taken into account. Let 
d denote the PPM type: d = 1 stands for the 3-PRR PPM; 
d^2 stands for the 3-RPR PPM; and <i = 3 stands for the 3- 
RRR PPM. 

As a result, the optimization problem contains one discrete 
variable, i.e., d, and five continuous design variables, i.e., R, r. 
Lb, rj and rp. Hence, the design variables vector x is given by: 

x^[dRrLhrjrpY (16) 

Formulation of the Optimization Problem 

The Multiobjective Design Optimization Problem of PPMs 
can be stated as: Find the optimum design variables x of PPMs 
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in order to minimize the mass of the mechanism in motion and 
to maximize its regular shaped workspace subject to geometric, 
kinematic and accuracy constraints. 

Mathematically, the problem can be written as: 



design variables and a discretization of this workspace is per- 
formed. The constraints of the optimization problem are also 
evaluated at each grid point of the regular-shaped workspace to 
check whether they are satisfied or not. A multiobjective genetic 



minimize f\ (x) = mx 
maximize f2{x)=R„ 

over X — [d R r Lt, rj rpY 



R 

subject to : gi:Lh + r>- 

§2 : < p; < V3R 
g3:K-'{J)>0.l 



(17) 



84 ■ k^y > 



F- 



oz 



: kT" > 



10 



8<j): 7r/180 

where x/^ and x„fo are the lower and upper bounds of x, respec- 
tively. 



RESULTS AND DISCUSSIONS 

The multiobjective optimization problem ([TT] ) is solved by 
means of modeFRONTIER 1241 and by using its built-in mul- 
tiobjective optimization algorithms. MATLAB code is incor- 
porated in order to analyze the system and to get the numer- 
ical values for the objective functions and constraints that are 
analyzed in modeFRONTIER for their optimality and feasibil- 
ity. The lower and upper bounds of the design variables are 
given in Tab. [T] The components of the PPMs are supposed 
to be made up of steel, of material density d = 7850 kg/m^ 
and Young modulus £ = 210x10'' N/m^. For each iteration. 



TABLE 1. LOWER AND UPPER BOUNDS OF THE DESIGN 
VARIABLES 



Design Variable 


d 


R[m] 


r[m] 


Lfc [m] 


i-j [m] 


'> [m] 


Lower Bound 


1 


0.5 


0.5 


0.5 








Upper Bound 


3 


4 


4 


4 


0.1 


0.1 



the regular-shaped workspace is evaluated for the corresponding 



TABLE 2. modeFRONTIER ALGOTITHM PARAMETERS 



Scheduler 


MOGA-II 


Number of iterations 


200 


Directional cross-over probability 


0.5 


Selection probability 


0.05 


Mutation probability 


0.1 


DNA (DeoxyriboNucleic Acid) string 


0.05 


mutation ratio 




DOE algorithm 


Sobol 


DOE number of designs 


30 


Total number of iterations 


30 X 200 = 6000 



algorithm (MOGA) is used to solve MOOP {Tji and to obtain the 
Pareto frontier in the plane defined by the mechanism mass and 
the workspace radius. modeFRONTIER scheduler and Design Of 
Experiments (DOE) parameters are given in Tab.|2] MATLAB is 
used to evaluate each individual of the current population (gen- 
erated by the modeFRONTIER scheduler). MATLAB returns the 
output variables that are analyzed by modeFRONTIER for the 
feasible solutions according to the given constraints. At the end, 
the Pareto-optimal solutions are obtained from the generated fea- 
sible solutions. 




200 400 600 800 1000 1200 1400 1600 
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FIGURE 6. PARETO FRONTIER OF MOOP 03 
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The Pareto frontier, solution of MOOP (fTTT l. is depicted in 
Fig.|6]whereas the design parameters and the corresponding ob- 
jective functions for two extreme and one intermediate Pareto 
optimal solutions, as shown in Fig. |6l are given in Tab. [3] The 
CAD designs illustrating the three foregoing solutions are also 
shown in Fig. [8] 

It appears that all Pareto-optimal solutions of MOOP ( fTTj i 
are 3-PRR PPMs. Accordingly, Fig. Q illustrates the Pareto 
Frontiers associated with the three planar parallel manipulator 
architectures. It is noteworthy that the Pareto-optimal solutions 
associated with the 3-PRR PPM architectures are better than 
the Pareto-optimal solutions associated with the 3-RPR and 3- 
RRR PPM ai-chitectures. 
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FIGURE 7. PARETO FRONTIERS ASSOCIATED WITH THE 3- 
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ARCHITECTURES 



Figures [9(a)| -(c) and |10(a)| -(c) show the evolution of the de- 
sign variables as a function of along the Pareto Frontier as- 
sociated with each PPM architecture. It is noteworthy that the 
higher R^., the higher the design variables. It is apparent that the 
variations in variables R, r, L/, and rj with respect to (w.rt.) R^. 
are almost linear whereas the variations in rp w.r.t. /?„, is rather 
quadratic. This is due to the fact that the higher the size of the 
mechanism the higher the bending of the moving platform links 
whereas the intermediate links are mainly subjected to tension 
and compression. 



CONCLUSIONS 

In this paper, the problem of dimensional synthesis of par- 
allel kinematics machines was addressed. A multiobjective de- 
sign optimization problem was formulated in order to determine 
optimum structural and geometric parameters of any parallel 
kinematics machine. The proposed approach is similar to that 



used in f25l but we took into account the mass and the regular 
workspace instead of considering the entire volume of the ma- 
nipulator The proposed approach was applied to the optimum 
design of three planar parallel manipulators with the aim to min- 
imize the mass in motion of the mechanism and to maximize 
its regular shaped workspace. Other performance indices can be 
used as constraints. However, they cannot necessarily be used 
as objective functions as the latter are usually formulated as a 
sum of an index over all the manipulator workspace. As another 
constraint, we could use the collisions between the legs of the 
manipulator 

REFERENCES 

[1] A. M. Hay and J. A. Snyman. Methodologies for the opti- 
mal design of parallel manipulators. International Journal 
for Numerical Methods in Engineering, 59(11):131-152, 
2004. 

[2] J. P. Merlet. Parallel Robots. Kluwer Academic PubUshers, 
Norwell, MA, USA, 2006. 

[3] Y. Lou, G. Liu, N. Chen, and Z. Li. Optimal design of paral- 
lel manipulators for maximum effective regular workspace. 
In Proceedings of the lEEE/RSJ International Conference 
on Intelligent Robots and Systems, pages 795-800, Alberta, 
2005. 

[4] Y. Lou, G. Liu, and Z. Li. Randomized optimal design of 
parallel manipulators. IEEE Transactions on Automation 
Science and Engineering, 5(2):223-233, 2008. 

[5] E. Ottaviano and M. relli. Workspace and optimal design 
of a pure translation parallel manipulator-tsai manipulator. 
Meccanica, 35(3):203-214, May 2000. 

[6] E. Ottaviano and M. Ceccarelli. Optimal design of capaman 
(cassino parallel manipulator) with prescribed workspace. 
In 2'"^ Workshop on Computational Kinematics KC2001, 
pages 35-43, Seoul, South Korea, 2001. 

[7] F. Hao and J. -P. Merlet. Multi-criteria optimal design of 
parallel manipulators based on interval analysis. Mecha- 
nism and Machine Theory, 40(2):157-171, 2005. 

[8] M. Ceccarelli, G. Carbone, and E. Ottaviano. Multi criteria 
optimum design of manipulators. In Bulletin of the Polish 
Academy of Sciences Technical Sciences, volume 53, 2005. 

[9] C. M. Gosselin and J. Angeles. The optimum kinematic 
design of a planar three-degree-of-freedom parallel manip- 
ulator ASME Journal of Mechanisms, Transmission and 
Automation in Design, 110:35-41, 1988. 
[10] C. M. Gosselin and J. Angeles. The optimum kinematic 
design of a spherical three-degree-of-freedom parallel ma- 
nipulator. Journal of Mechanisms, Transmissions and Au- 
tomation in Design, 1 ll(2):202-207, 1989. 
[11] H. H. Pham and I-M. Chen. Optimal synthesis for 
workspace and manipulability of parallel flexure mecha- 
nism. In Proceeding of the If'' World Congress in Mech- 



8 



Copyright © 2010 by ASME 



TABLE 3. THREE PARETO OPTIMAL SOLUTIONS 
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